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Abstract. In the theory of thermonuclear reaction rates, analytical evaluation of thermonuclear 
functions for non-resonant reactions, including cases with cut-off and depletion of the tail of the 
Maxwell-Boltzmann distribution function were considered in a series of papers by Mathai and 
Haubold (1988). In the present paper we study more general classes of thermonuclear functions 
by introducing a pathway parameter a, so that when a — > 1 the thermonuclear functions in 
the Maxwell-Boltzmannian case are recovered. We will also give interpretations for the pathway 
parameter a in the case of cut-off and in terms of moments. 

1. Thermonuclear Functions 

The standard thermonuclear function in the Maxwell-Boltzmann case in the theory of nuclear 
reactions, is given by the following (Critchfield, 1972; Haubold and Mathai, 1985; Mathai and 
Haubold, 1988): 



h= x^e-^-^dx, a > 0, b > 0, p > 0, x > 0. 
J o 
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Several modifications to this standard thermonuclear function in (1) are considered in Mathai and 
Haubold (1988), Saxena et al. (2004), and Lissia and Quarati (2005). One such case is the cut-off 
of the distribution function at point d, in which the integral becomes 

xJ -i e -ax-bx-" dx ^ a > o, 6 > 0, p > 0, d < oo. (2) 
Another case is depletion of the tail of the distribution function where the integral becomes 

oo 

xJ -i e -ax 5 -bx-" dx ^ a > o, 6 > 0, p > 0, (5 > 0. (3) 

Note that if S = 1 is taken as the standard Maxwell-Boltzmannian behavior, then for 8 > 1 the 
right tail will deplete faster and if 5 < 1 then the depletion will be slower in I3. We consider more 
general classes of (1), (2), and (3) by replacing e~ ax by a binomial factor [1 — a(l — a)x s ] 1 - a so 
that the integrand x 7_1 [1 — a(l — a)x 5 ] 1 ~ a stays in the generalized type-1 beta family when a < 1, 
type-2 beta family when a > 1 and generalized gamma family when a — > 1. Thus we consider the 
general class of integrals 



r(<5) 
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c 7_1 [l - a(l - a)x 6 



-bx 



dx, 
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for a > 0, 5 > 0, b > 0, p > and 



f d 

/W= / x 7 - 1 [l-a(l-a)x' 5 ]^e- 6 ^ P dx, 
./o 
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for a > 0, d > 0, 5 > 0, b > 0, p > 0. 

For a < 1 we have extensions of generalized type-1 beta family of integrals in (4) and (5), 
whereas for a > 1 we have extensions of generalized type-2 beta family of integrals in (4) and 
(5). When a —> 1 both of these forms will go to I 3 in (3). This consideration is motivated by 
particularly taking into account Tsallis statistics as developed in Tsallis (1988, 2004), Gell-Mann 
and Tsallis (2004), and discussed in Cohen (2005). 

The pathway parameter a enables one to go to different functional forms in (4) and (5), thus 
creating a pathway through these different families of functional forms as shown in Section 2. Also 
in Section 2 a connection to generalized entropy of order a is given. In Section 3 interpretations for 
a in terms of the cut-off d in (2) as well as in terms of moments are given. Section 4 evaluates the 
thermonuclear functions in terms of H-functions. Moments are evaluated in Section 5. Conclusions 
are drawn in Section 6. For all Sections, concerning generalized entropies of order a, see Mathai 
and Rathie (1975), Gell-Mann and Tsallis (2004), and Mathai and Haubold (2007). 

2. Optimization of an Entropy 

Let us start with the optimization of the generalized entropy of order a, 

Jx « - 1 

under the conditions 

(i): f(x) > for all x, 
(h): f x f(x)dx < oo, 

(iii) : j' x x^ l - a ^f{x) dx = fixed for all /, 

(iv) : J x xt( 1_q )+' 5 f(x)dx = fixed for all /, where 7 and 8 are fixed parameters. 

In this discussion we will consider 1 to be a real scalar variable. For fixed a, optimization 
of M a (f) implies optimization of J x [f(x)] 2 ~ a dx subject to conditions (i) to (iv). If calculus of 
variation is used then the Euler equation is the following: 

^j[f~ a - Xrx^-^f + X 2 x^ 1 -^+ 5 f] = => 

(2 - a)/ 1 "" = X 1 x^ 1 -^ [1 - ^V], ajtl,2,=> 

f x = ClX i[i - a (i - a )x 5 ] ^ (7) 

for ^ = a(l — a) for some a > 0, c = (2"^) 1-01 • Here c\ can act as the normalizing constant when 
/ in (7) is a statistical density function. 

For a<l,a>0,(5>0in(7)we need an additional condition to make (7) a density, namely 
1 — ail — a)x s > < x < — r when x is positive. When a > 1 then writing 1 — a = 

[o(l-a)]7 

— (a — I) we have 

f 2 = c 2 x' r [l + a(a-l)x 5 }-^,a> l,a>0,<5>0. (8) 
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When a — > 1, both fi and /2 go to /3, where 

/ 3 = c 3 rc 7 e- Q ^ , a > 0, (5 > 0. (9) 



„<5 



In (7), (8) and (9) one can interpret 8 as the rate of depletion of the right tail. If 8 = 1 is taken 
as the standard behavior, then for 8 > 1 the tail is depleted at a faster rate and when 8 < 1 it is 
depleted at a slower rate. Then 8 measures the rate of depletion of the right tail. Note that the 
normalizing constants ci, C2, C3 will be different for the three cases — 00 < a < l,a > 1 and a < 1. 

The scalar version of the pathway model of Mathai (2005) is 

/ = c \xy[l - a(l - a)\x\ 5 ]^,v > 0,a > (10) 

for — 00 < x < 00. Our cases in (7), (8), (9) are special cases of (10) for 77 = l,x > 0. Note that 
the normalizing constant c in (10) is also different for the three cases a < l,a > 1, and a — > 1. 
Observe also that in (8) and (9), 5 could be negative also. In (10) for a > 1 and a — > 1, 5 could 
be negative also. 

We look at the normalizing constants ci, 02,03 in (7) to (9) for x > 0. For a < 1 make the 

1 _ 1 

substitution y = ail — a)x 5 =. Then dx = — — — -^t- and then 

5[o(l-a)]7 

/i(x)dx 

- / y^ -1 (i - y) 1 ^?/ 
■/ 



5[a(l-a)]- 



rp^rfe + 1) 



"*[a(l-a)] 2 * 1 r( T ^ + 1 + 2+1) 
evaluating the integral with the help of a type-1 beta integral, for 7 + 1 > 0, a < 1. Then 

7+1 r( T 1 - + i + 4i) 

Cl = 8[a(l - a)} 2 ^ : 5 — (11), 

r(2±i)r(^- + i) 

for a < 1, 7 + 1 > 0, 8 > 0, a > 0. The h-th moment of x, denoted by E(x h ), where E denotes the 
expected value, is then given by 

( ] [a(i-a)]* r ( ^±i) r (T ^ + i + ^)' ( j 

for a<l,7 + l + /i>0, <5 > 0, a>0. Note that all the h-th moments for h > — 1 exist when 7 = 0. 
Proceeding the same way we have the following h-th moment for a > 1. That is, for a > 1, 

E{x h ) = ( ' ] U -j ^-i (13) 

[a(a-l)]* r(2+l) r(^-i±i) 

for ^-2±1 >0,^-2± 1 ±ii >0, 7 + l + /i> 0,7 + 1 >0 or -7 - 1 < h < ^ - 7 - 1, a > 1, 
and the normalizing constant 

c 2 = < y[ a (a- l)] 2 * 1 r( "T l) Ti — ( 14 ) 

r(2±i)r(^-^±i) 
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for ^-j- — :! ^- >0,7 + l>0, 5 > 0, a > 1. When a — > 1 then c\ and C2 reduce to 

C3 = f(Wy' 7 + 1>0 ( 15 ) 



and for a — > 1, 



E(x h ) = ± (16) 



for a— >l,7 + l + /i>0,7 + l>0, a > 0, <5>0. Proceeding the same way one can compute c and 
E(\x\ h ) in (10) for the three cases a < 1, a > 1 and a — > 1. Since the procedure is parallel we will 
not give the details here. 

3. Pathway Parameter a 

In the above formulation the parameter a creates a pathway to three different functional forms. 
For a < 1 the pathway density of (10) stays in the generalized type-1 beta family of distributions. 
As a moves closer and closer to 1 we move continuously to a generalized gamma or exponential 
type family of densities. As a goes above 1 we move from the generalized gamma type family 
to a generalized type-2 beta form. Thus the pathway parameter a takes one to three different 
functional forms. This is the distributional pathway. 

Let us see what happens to the generalized entropy of order a in (6), when a moves from — oo 
to oo. When a approaches 1, M a (f) comes closer to Shannon's entropy and finally 

lim M a (f) = — f(x)\nf(x)dx = S(f) = Shannon entropy. 

a->l J 

Then when a increases from 1, M a (f) moves away from Shannon entropy to the form in (6). Thus 
a creates a pathway for the generalized entropy in (6) also. This is the entropic pathway. 

The differential pathway provided by the pathway parameter a can be seen from the following: 
As an example, let us start from (2). Let 

g (x) = = z 7 [l - a(l - a)x 5 }^,x > 

Cl 

± g ( x ) = lg( x )-aSx 5 - 1+ ^[g(x)r 

= -a[g(x)r for 7 = 0,5=1. (17) 

This is the power law. The differential equation satisfied by gs = ^ of (9) is given by 



^-g 3 (x) = ( 1 -a5x 5 - 1 )g 3 (x) 
ax x 

= -ag 3 (x) for 7 = 0,5 = 1. (18) 

Thus when a moves to 1 the differential pathway is from the power law in (17) to the Maxwell- 
Boltzmann law in (18). 

3.1. An interpretation of the pathway parameter a 
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Let us start with the pathway density in (7). For this to remain a density we need the condition 
1 — ail — a)x 5 > or when x is positive then < x < — r , a < 1, or if we have the model in 

[a(l-a)]7 

(5) then 

1 1 

r < x < r , a < 1. 

[a(l-a)]i [a(l-a)]* 

Outside this range the density is zero. Thus for x > the right tail of the distribution function is 
cut-off at — r . If d is this cut-off on the right then 

[o(l-a)]7 

i 1 

a = 1 -tt for a < 1. (19) 



[a(l-a)]* ' ad 5 

As a moves closer to 1 then the cut-off d moves farther out and eventually when a — > 1 then 
c£ — > oo. In this case a is computed easily from the cut-off. 

4. Evaluation of the Generalized Thermonuclear Functions 

Evaluation of l[ 5 J and I^J in (4) and (5) can be done by using the convolution property of 
Mellin transforms by writing in the form 

[°°-fi(x)M-)dx 

J X X 

where 

/i(x) = x 7 [l - a(l - a)x 5 ]^ and / 2 (x) = e _a;P (20) 

so that /2(f) = e~" Pa; P and that b~p = u. For a > 1 we write /i(x) = x 7 [l + a(a — l)x <5 ]~ 
and for a < 1 we write /i(x) = x 7 [l — o(l — a)x ] l - a . According to the convolution property of 
Mellin transforms, the Mellin transform of is t ne product of the Mellin transforms of f\ and 
/2, denoted by Mf 1 (s) and Mf 2 (s) respectively, where for a > 1, taking the more generalized forms 
of /1, we have for 5 > 

/>oo 

M fl (s)= x 7+s " 1 [l + a(a- l)x 5 ]-^rdx 

Ka-l)]^^^) «-l ^ 



/•OO 1 „ 

M f Js)= x s - 1 e~ xP dx= -r(-),s>0. 



Thus for a > 1, 



for s > 0, 7 + s > 0,p > 0, 5 > 0, a > 0, 
available by taking the inverse Mellin transform. That is, 



1 ^T 1 > °' a > L Then the integral in l[ 5 J is 



AS) 

'la 



Ma(«-i)]lr(^) 



— / r(-)r(- 

2 ™ J c - ioo p 



-)r( 



(7 + Sy 



a- 1 



)(ftp[a(a- l)]*)" s ds 



1 -tf?' 1 



a r / i \ "1,2 



Ma( a -i)]ir(^) 



Ha^-l)]*^^ 



, -y In' 

i \ (a-1) 



1 u-2,0 



l I 

a 5 op 



, for a — > 1 



(21) 



For the definition, theory and applications of H- functions see Mathai and Saxena (1978). Note that 
in Mf 1 (s) we could have considered the form with x~ 5 , 5 > instead of x 5 , 5 > 0. Mathematically, 
both integrals are tractable. For x s with 5 = 1, we have the extension of the Maxwell-Boltzmann 



case e 



For a < 1 we may evaluate the pathway integral corresponding to ijjfj , namely, 



l£ = [ x^[l-a{l-a)x 5 ]^ 
J o 



e bx P dx, a < 1 



where for a < 1, 



fi (x) = x 1 [1 — a{l — a)x S ] 1 - a dx 



for d - 
Then 



- — t or < x < — r and fi(x) = elsewhere, and /2(f) = e xP ,0 < x < 00. 



[o(l-a)]7 



[a(l-a)]7 



rd 1 

M /l (a)= / x 7+s_1 [l -a(l -a):/]~, d= -,a<l 

[a(l — a)]' 



r(^)r(^ + i) 



5[a(l-a)] I ^r( T ^ + l + ^) 

/•oo Is 

M h (s)= / x^e^'ds = -T(-),8 > 0,p > 0. 
"Jo P P 



,7 + s>0,a<l 



Hence, 



/-(*) 
'2a 



1^ 1 1 r +io ° r ( l2 T £l ) r (T^ + i )„i 



6 P T{ P ] [a(l - a)]i 2iri J c _ loo r( _l_ + 1 + (l+£i) 



(&p[a(l -a)]^)- s ds 



p£[a(l - a)]~ 

L_ w 2 >° 

"t "0,2 



7 "1,2 



bp [a(l — a 



7 1 1 

op a 5 



(o,i),(*,i) 
, for a — > 1. 



For 5 = 1 we have the extension of the Maxwell-Boltzmann case of e aa: in fi(x). 
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5. An Evaluation of a 



In this connection let us evaluate the h-th. moment in the density 

ft(x) = cix 7 [l + a(a - l)x 5 ]"^ T ,a > 1,6 > 0,a > 0,x > 

where c\ is the normalizing constant. Integration gives the following, denoting the h-th. moment 
by E(x h ). 

i r(&m) r(^ - Mjtil) 



£(x h ) 



[a(a-l)]* r(^±i)) r(^-(^) ' 



for a > 1,6 > 0,j + h + 1 > 0, ^ - (7+ ^ +1) > 0, 7 + 1 > 0. When § = m,m = 1,2,... or when 
/i = 6m, m = 1, ... we have for m = 1, 



which gives 



1 (7 + 1) 



a(a-l) 5 i x 

a— 1 <5 

( 7 + 1 + 25) ( 7 + 1) 



(7+1 + 5) (7 + 1 + 6)aE(x s ) 
for 1 < a < 1 + . Similarly one can evaluate a for a < 1 also from E(x s ). In this case 

1 (7 + 1)! 

a = 1 + 7~~ [o 



(5 + 7 + 1) L a£(x 5 ) J 

for 7 + 1 > a6E(x 5 ). 
6. Conclusions 

In the field of stellar, cosmological, and controlled fusion, for example, the core of the Sun is 
considered as the gravitationally stabilized solar fusion reactor. The probability for a thermonuclear 
reaction to occur in the solar fusion plasma depends mainly on two factors. One of them is the 
velocity distribution of the particles in the plasma and is usually given by the Maxwell-Boltzmann 
distribution of Boltzmann-Gibbs statistical mechanics. The other factor is the particle reaction 
cross-section that contains the dominating quantum mechanical tunneling probability through a 
Coulomb barrier, called Gamow factor. Particle reactions in the hot solar fusion plasma will occur 
near energies where the product of velocity distribution and tunneling probability is a maximum. 
The product of velocity distribution function and penetration factor is producing the Gamow peak. 
Mathematically, the Gamow peak is a thermonuclear function. In case of taking into consideration 
electron screening of reactions in the hot fusion plasma, the Coulomb potential may change to 
a Yukawadike potential. Taking into account correlations and long-range forces in the plasma, 
the Maxwell-Boltzmann distribution may show deviations covered by the distribution predicted by 
Tsallis statistics in terms of cut-off or depletion of the high- velocity tail of the distribution function. 
In this paper, closed-form representations have been derived for thermonuclear functions, thus for 
the Gamow peak, for Boltzmann-Gibbs and Tsallis statistics. For this purpose, a generalized 
entropy of order a and the distribution function, emanating from optimizing this entropy, has 
been considered. The case a = 1 recovers the Maxwell-Boltzmannian case. This general case 
is characterized by moving cut-off, respectively the upper integration limit of the thermonuclear 
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function to infinity. The closed-form representations of thermonuclear functions are achieved by 
using generalized hypergeometric functions of Fox or H-functions. 
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